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Abstract 


The  problem  of  detecting  a  signal  known  except  for  amplitude  in  non- 
Gaussian  noise  is  addressed.  The  noise  samples  are  assumed  to  be 
independent  and  identically  distributed  with  a  probability  density  function 
known  except  for  a  few  parameters.  Using  a  generalized  likelihood  ratio 
test  it  is  proven  that  for  a  symmetric  noise  probability  density  function 
the  detection  performance  is  asymptotically  equivalent  to  that  obtained  for 
a  detector  designed  with  a  priori  knowledge  of  the  noise  parameters.  A 
computationally  more  efficient  but  equivalent  test  is  proposed  and  a 
computer  simulation  performed  to  illustrate  the  theory. 
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I,  Introduction 


A  large  body  of  knowledge  is  available  in  statistical  detection 
theory  for  the  detection  of  a  known  deterministic  signal  in  noise  with 
a  known  probability  density  function  (PDF) * ^  In  practice,  the  noise 
characteristics  are  never  known  a  priori  since  they  depend  on  unknown 
or  incompletely  understood  physical  phenomenon . ^  Relatively  little 
is  known  about  the  design  of  optimal  detectors  for  this  situation.  One 
approach  which  has  been  investigated  is  to  use  Bayesian  methods  and 
assign  priors  to  the  unknown  parameters  of  a  noise  PDF.  L  1  The 
resulting  "optimal*  detector  requires  multidimensional  integration 
which  is  generally  not  practical.  Furthermore,  the  performance  will  be 
critically  dependent  on  the  choice  of  the  priors;  the  appropriate 
choice  is  never  known  a  priori.  The  entire  problem  becomes  much  more 
difficult  when  the  noise  can  no  longer  be  characterized  as 
Gaussian. t81.t91.110] 

The  approach  which  is  taken  here  is  to  apply  the  theory  of 
generalized  likelihood  ratio  testing  for  composite  hypothesis 
testing.  ^  The  noise  PDF  is  assumed  to  be  known  except  for  a  finite 
set  of  parameters.  The  unknown  parameters  are  then  estimated  using 
maximum  likelihood  estimators.  This  approach  does  not  require  the 
arbitrary  selection  of  priors  as  in  the  Bayesian  detector.  The 
difficult  multidimensional  integration  required  by  the  Bayesian 
detector  is  replaced  by  a  multidimensional  maximization  of  the 
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likelihood  function.  Of  course,  this  maximization  may  be  also 
difficult  to  implement.  By  retaining  the  parametric  form  of  the  noise 
PDF  it  is  expected  that  the  detection  performance  will  exceed  that  for 
the  optimal  nonparametric  or  robust  detector. 

Accepting  the  desirability  of  a  generalized  likelihood  test 
(GLRT),  several  important  properties  are  shown  to  hold  for  large  data 
records.  The  main  theorem  states  that  the  asymptotic  (large  data 
record)  performance  of  the  GLRT  is  equal  to  that  of  a  clairvoyant  GLRT 

detector  if  the  noise  PDF  is  symmetric.  (Examples  of  symmetric  PDFs 
are  the  Gaussian,  Laplacian  and  Gaussian  mixtures.)  The  claivoyant 
GLRT  is  one  which  designs  the  detector  using  perfect  knowledge  of  the 
unknown  noise  parameters.  Hence,  the  GLRT  asymptotically  achieves  an 
upper  bound  in  performance  and  so  can  be  said  to  be  optimal .  The 
theoretical  results  are  quite  general  and  apply  to  many  practical  cases 
of  interest.  Computer  simulation  results  confirm  the  asymptotic  theory 
for  finite  data  records.  Similar  but  less  general  results  had  been 
obtained  by  the  author  in  [12]. 

The  paper  is  organized  as  follows.  Section  II  provides  a  review 
of  the  GLRT  and  its  asymptotic  properties.  The  GLRT  as  applied  to  the 
detection  problem  is  described  in  Section  III  while  its  performance  and 
the  main  theorems  are  given  in  Section  IV.  Section  V  applies  the 
theory  to  detection  in  non-Gauss ian  noise  while  Section  VI  describes 
the  resu’ts  of  a  computer  simulation. 
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II  Review  of  Generalized  Likelihood  Ratio  Testing 


This  section  summarizes  the  important  theory  to  be  used  later. 
The  material  has  been  extracted  from  references  [1] , [2] , [3] .  Consider 
the  problem  of  testing  the  value  of  a  k  =  r+s  dimensional  parameter 

e 

-r 

e 

“S. 

where  0r  is  rxl  and  0g  is  sxl.  A  common  hypotheses  test  is  that  9  lies 
in  an  r  dimensional  subspace  or 


H0: 

— r 

-  ?•  «s 

(1) 

Hi: 

ir 

*  2-  «, 

0g  is  of  no  concern  and  may  take  on  any  value.  It  is  sometimes 
referred  to  as  a  nuisance  parameter.  Assuming  the  data 
x  =  [xq  x«^  xN-l]  are  °^serve^  with  a  joint  probability  density 

function  p ( x  ;  0r,  0S ) ,  a  generalized  likelihood  ratio  test  for  testing 
(1)  is  to  decide  if 


A  A 

p(x  ;  0  ,9  ) 

L  =  ~r  ~s 

A 

p(x  ;  0,§s) 


>  r 


(2) 


A 

for  some  threshold  y.  9§  is  the  maximum  likelihood  estimator  (MLE)  of 

A  ft 

0g  assuming  Hq  is  true  while  0r,  0g  are  the  joint  MLEs  of  0r  and  0g 

A 

assuming  is  true,  0$  is  found  by  maximizing  p  (x  ;  0,9S)  over  0g, 
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Likewise,  0  ,  0  are  found  by  maximizing  p(x  j  0„,  0O)  over  0  .  0„. 

The  statistics  of  Lq  are  difficult  to  obtain  in  general.  For 
large  data  records  (N  large)  or  asymptotically  it  may  be  shown  that 
2  IiiLq  is  distributed  in  the  following  manner, 


2  InL  ~  X  2  under  HA 

b  r  0 

2 

2  InL  ~  X'  ( X,r )  under  H. 

b  1 


(3a) 

(3b) 


Here  Xr2  represents  a  Chi  squared  distribution  with  r  degrees  of 
freedom  and  X'2(X,r)  represents  a  noncentral  Chi  square  distribution 
with  r  degrees  of  freedom  and  noncentrality  parameter  X,  Note  that 
X'  (0,r)  —  Xr2  or  the  distribution  under  Hq  is  a  special  case  of  the 
distribution  under  and  occurs  when  X  =  0.  The  noncentrality 

parameter  X,  which  is  a  measure  of  the  discrimination  between  the  two 
hypotheses,  is 


x  ■  £  [j«rer(2-2,)  -  I,res <»•«»>  I^es<5-2*,4e,<?-?s)]  2' 


(4) 


where  0r,  0g  are  the  true  values.  The  terms  in  the  brackets  of  (4)  are 
found  by  partitioning  the  Fisher  information  matrix  for  0 


1(0)  = 


-0r0r 

(0r, 

0S) 

V.‘Sr- 

es)  ' 

~0S0r 

(0r. 

0S) 

(§r* 

0S) 

and  the  partitions  are  defined  as 


(5) 
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r  x  r 


All  the  partitions  of  the  Fisher  information  matrix  are  evaluated  at 
0^  =  0  and  the  true  value  of  for  use  in  (4), 

The  motivation  for  using  a  GLRT  is  that  for  large  data  records  it 
exhibits  certain  optimality  properties.  (It  should  be  noted  that  in 
general  the  hypothesis  testing  problem  of  (1)  does  not  admit  a 
uniformly  most  powerful  (UMP)  test.^)  These  optimality  properties 
are  that  of  all  the  tests  which  are  invariant  to  a  natural  set  of 
transformations  the  GLRT  exhibits  the  largest  probability  of  detection 
for  all  values  of  the  unknown  paramaters.  The  GLRT  is  said  to  be  the 
asymptotically  uniformly  most  powerful  invariant  (UMPI)  test.^  The 
conditions  under  which  the  asymptotic  results  apply  to  finite  length 
data  records  are  difficult  to  quantify  in  general.  Heur ist ically ,  it 
may  be  said  that  the  asymptotic  results  will  apply  if 

1)  The  MLEs  necessary  to  define  the  GLRT  are  adequately  characterized 
by  their  asymptotic  properties  or  the  PDF  of  the  MLE  is  Gaussian 
with  mean  equal  to  the  true  parameter  value  and  covariance  matrix 
equal  to  the  inverse  of  the  Fisher  information  matrix  and 
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2) 


The  value  of  0r  when  is  true  is  close  to  the  value  when  Hq  is 
true  or  0r  ~  0.  In  effect,  it  says  that  we  must  be  testing  for 
slight  departures  of  0r  from  zero. 

To  implement  the  GLRT  requires  one  to  find  the  MLEs  under  Hq  and 
In  many  cases  this  is  a  difficult  analytical  task.  To  avoid  some 
of  the  difficulties  use  can  be  made  of  the  Rao  test  which  is 
asymptotically  equivalent  to  the  GLRT.  The  Rao  test  ^  decides  Hj 

if 


V 


dlnp(x  ;  0  ,0  ) 
—  — r  — s 


ae 

— r 


A  31np (x  ;  0  ,0  ) 

V  (0,9  )  - ^ - - — — 

6  =  0  ‘  "  "s  r 

— r  - 


>  n  (7) 


e  =  o 

-r  - 


where 


V(0  ,  9  )  = 
-  -r  ~s 


h  9  -  h  9  ^9  9<§r-§s> 


r  r 


r  s 


s  s 


r  s 


Note  that  only  the  MLE  of  0g  under  Hq  need  be  found.  The  MLE  of  §r  and 
0$  under  H^  are  no  longer  required.  Also,  the  statistics  of  LR  are 
given  by  (3  ) . 


Ill  Statement  of  Problem  and  Detector 


We  will  consider  the  problem  of  detecting  a  signal  known  except 
for  amplitude  in  independent  and  identically  distributed  non-Gaussian 
noise.  The  univariate  PDF  of  the  noise  is  assumed  to  be  known  except 
for  a  finite  set  of  parameters.  Mathematically,  we  have 
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(8) 


H0:  xt  =  nt 

H1:  Xt  =  ^St  +  nt 
for  t  =  0,  1,  . N  -  1.  p,  the  amplitude,  is  unknown  and  can  take  on 

any  value  while  st  is  known.  The  PDF  of  nt,  which  is  the  same  for  any 

t  is  denoted  by  p(n;  a),  where  a  is  a  set  of  unknown  parameters.  The 

PDF  for  the  noise  is  not  restricted  to  be  Gaussian,  only  that  it 

satisfies  certain  regularity  conditions  to  insure  the  validity  of  the 

rq  ] 

asymptotic  MLE  properties. 

The  problem  of  (8)  can  be  recast  as  a  test  of  the  amplitude 


Hq:  p.  =  0,  a  (9) 

H^:  p  7^  0,  a 


which  is  recognized  as  a  composite  hypothesis  testing  problem  of  the 
form  of  (1).  Because  p  may  take  on  positive  or  negative  values  it  is 
well  known  that  no  UMP  test  exists,  even  for  a  priori  knowledge  of 
a.^  Hence  a  GLRT  is  proposed.  Due  to  the  independence  of  the  noise 
the  GLRT  of  (2)  reduces  to 


LG 


T  p(x  “  ps  ;  a) 
t=0  z  Z 


N-l 

T  p(x  ;  a) 
t=0 


(10) 


where  the  identification  0  =  [  0rT  0*  ]T  =  [  p  aT  ]T  has  been  made. 


p  is  the  BILE  of  p  under  and  a,  a  are  the  MLEs  of  a  under  and  Hq, 


respectively.  p ,  a  are  found  by  maximizing 
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N-l 

T  p(x  -  ns  ;  a) 

t=0 

A 

over  | x  and  a  while  a  is  found  by  maximizing 
N-l 

T  p(x  ;  a) 
t=0  r 

over  a. 


IV  Asymptotic  Performance  of  GLRT 


The  asymptotic  performance  of  the  GLRT  given  in  (10)  may  be  found 
by  making  use  of  the  results  of  Section  II.  Making  the  identifications 


so  that  0r  is  a  scalar  or  r  =  1  and  a  is  an  sxl  vector,  we  have  from 
(3), (4) 


2  IiiLq  -  Xjl*  under  Hq  (11) 

-  X'2  (1,  X)  under 


where  X  =  n*  [  I  (0,  a)  -  I  (0,  a)  I  X(a,  a)  IT  (0,  a)  1  (12) 

L  “  -pa  -  -aa  -  -  -pa  -  J 
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The  probability  of  detection 


PD  =  Pr  {  21nLG  >  y'  1^} 

and  the  probability  of  false  alarm 


(13) 


PFA  “  Pr  {  21nLG  >  T*  |  H0J 


(14) 


are  easily  found  by  noting  that  a  X'2(1,X)  random  variable  is 
equivalent  to  the  square  of  a  normal  random  variable  with  mean  and 
variance  1.  It  can  be  shown  that 


pfa  = 


(15) 


pd  -  qc^y'  -  ^1)  +  q  (^  +  >11) 


(16) 


1  2 

As  1  i  ^ 

where  Q  (£)  53  L  —  -  ■■  e  dx 

\|2 n 

It  is  interesting  to  note  that  the  probability  of  detection  increases 
monotonically  with  X.  Furthermore  a  special  case  of  the  theory 
presented  in  Section  II  occurs  when  0  ,  the  nuisance  parameters,  are 
known.  For  the  problem  at  hand  0g  =  a.  If  a  is  known,  then  the 
clairvoyant  GLRT  as  given  by  (10)  becomes 
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L, 


N-l 

T 

t=0 


(17) 


GC 


N-l. 

T  p(x  j  a) 
t=0 


and  the  asymptotic  statistics  of  2  lnLgp  are  given  by  (11)  but  with 


(18) 


Note  that  the  performance  of  the  GRLT  when  a  is  known  is  always  better 
than  or  at  worst  equal  to  the  performance  when  a  is  unknown  since 


Xc  1  X.  The  loss  in  performance  when  a  is  unknown  may  be  attributed  to 


the  reduction  in  Fisher  information  due  to  having  to  jointly  estimate  jjl 
and  a,  even  though  a  is  of  no  concern.  If  it  happens  that  =  X  then 

the  GLRT  will  produce  the  same  detection  performance  as  the  clairvoyant 

GLRT  detector,  i.e«,  the  GLRT  which  incorporates  a  priori  knowledge  of 

the  unknown  parameter  a.  From  examination  of  (12)  and  (18)  this  will 
be  the  case  if 


(19) 


This  condition  is  equivalent  to  requiring  the  MLEs  of  ji  and  a  to  be 
asymptotically  uncorrelated  (and  hence  independent)  when  ji  is  close  to 
zero . 

An  example  which  illustrates  this  behavior  is  the  following 
classical  problem.  Consider  the  detection  problem  of  (8)  with  st  =  1 
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and  ~  N(0,  cr2 ) .  The  noise  samples  are  independent  and  identically 
distributed  and  the  variance  a2  is  unknown.  The  amplitude  ji  can  be 
either  positive  or  negative  so  that  no  UMP  test  exists.  The  nuisance 
parameter  in  this  problem  is  the  noise  variance  a2.  If  one  derives  the 
GLRT,  (see  Appendix  A)  then  the  test  statistic  is 


2 InL  =  Nln  1  - 
(j 


-N-l 
1  v  2 

NJo  Z* 


1 N_1 

where  x  =  rT  I  x  ^ 
N  t=o 


which  is  the  sample  mean. 


(20) 


If  however  a2  were  known,  then  the  GLRT  would  yield 


21nLGC 


2 


X 


(21) 


which  might  be  referred  to  as  a  clairvoyiant  detector  since  it  makes 

2 

use  of  a  priori  knowledge  of  or  .  Now  assume  N  — > 00  so  that  by  the  law 
of  large  numbers  x  -> p  and  assume  n  is  close  to  zero.  Then  from  (20) 
21nL,Q  converges  to 


21nL(3 


Nx 


,  N-l 
1  T  * 

N  ^  Xt 
t=0  C 


(22) 


using  the  relation  ln(l-x)  ~  x  for  x  <<  1  .  But  by  Slutsky's  Theorem 
2lnLg  converges  in  distribution  to  that  of  21uLq£.  Recall  that  it  is 
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assumed  jx  — >0  and  hence 


N-l 

t=0 


2 


2 

X 

t 


a 


2 


whether  or  not  a  signal  is  present.  Since  the  distributions  of  Lq  and 
Lgc  are  asymptotically  identical,  their  performances  are  identical.  It 
is  an  easy  matter  to  show  that  (19)  holds  or 

Ilia2(  °>  °2)  =  0  (23 ) 

(Actually  (|±*  cr3 )  =  0  for  any  |i  in  this  case.)  A  careful 

examination  of  the  form  of  I^a*(p,  ff2 )  reveals  that  the  symmetric 
nature  of  the  Gaussian  PDF  is  responsible  for  (23).  In  fact  as  is  now 
proven  I  (  \l,  a)  will  be  zero  whenever  the  PDF  of  the  noise  is  an  even 

r  ^ 

function. 


Theorem  1:  Assume  that  N  samples  of  a  process  x^  =  \is ^  +  n^  are 

observed,  where  the  noise  samples  are  independent  and  each  one  is 

distributed  according  to  the  PDF  p(n;  a)  which  depends  on  the 

parameter  a.  s^,  the  signal,  is  known.  If  p(-n;  a)  =  p(n;  a),  for  all 

a,  then  I^a(^i,  a)  =  2'  where  the  Fisher  information  is  based  on 

{ x0  *  xl '  • • • *  XN-1  ^ • 
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Proof:  By  the  definition  of  the  Fisher  information 


I  (p,  a) 
-pa  - 


E 


N-l 

d In  T  p(x  ”  ps  ;  a) 
t=0 

d\i 


N-l 

din  T  p(x  -  ps  ;  a) 
t=0 

da 


(24) 


Since  the  noise  samples  are  independent,  p(x^  -  ps^;  5^  is  independent 
of  p(xj  “  psjj  ?)  for  j  ^  i  and  each  term 


dlnp (x^  -  ps^; 

30 

where  0  =  p  or  a^  has  expectation  zero.  Hence, 


I  (p,  a) 
-pa  r  - 


N-l 

1 

t=0 


91np(xt  -  ps^;  a) 

31np(xt  -  pst;  a)  1 

9p 

da  J 

(25) 


Consider  a  single  term  in  the  summation,  which  may  be  written  as 

f 


I- 


dlnp(xt~  \is^;  a)  dlnptx^-  ps^;  a) 


3p 


da 


p(xt  -  pst;  a)dxt 


Using  the  chain  rule  this  becomes 


-s , 


00  dlnp(x^-  p s^;  a)  dlnp(x^-  jis^;  a) 


dUt  -  |ist) 


da 


p(  x  -  (as ^ ;  a)dxt 


Letting  =  xt  -  pst,  we  have 
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00  dlnpCn^;  a) 


31np(nt;  a) 
0a 


p(n^;  a)dnt 


That  this  last  expression  equals  zero  follows  from  the  assumption  that 
p(nt;  a)  is  an  even  function  in  n^..  For  p(nt;  a)  even 


dlnpCn^;  a) 
0a 

is  even  in  nt  and 


01np (n 
0n 


1 1 

t 


a) 


is  an  odd  function  in  nt*  The  integral  is  then  an  odd  function  which 
is  integrated  over  an  even  interval  resulting  in  zero*  To  prove  the 
even  property  of  01np(n^;  a)/0a  consider  the  ith  component 


01np  (n^_ ;  a) 

0a . 
i 


lim 
A  0 


lnp(n^;  a  +  A^)  -  lnpCn^s  a) 

T. 

i 


i 

where  A^  is  a  vector  with  all  zeros  except  for  the  ith  element  which  is 
the  small  increment  A^.  It  is  clear  from  examination  of  the  partial 
derivative  definition  that  if  the  PDF  is  an  even  function  for  all  a, 
then  01np(nt;  a)/0a^  is  also  even.  Next  consider 
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dlnpCn^;  a) 
dn 


31np(-nt;  a) 
3n 


1  im 
A  — >  0 

1  im 
A  — >0 

1  im 
A 


lnp(n^  +  A;  a)  -  lnp(n^;  a) 


lnp(-nt  +  A;  a)  -  lnp(-n^;  a) 


lnp (n^  -  A;  a)  -  lnptn^;  a) 


-  1  im 
A  0 


lnpCn^;  a)  -  lnp(n^  -  A;  a) 


dlnp(n^;  a) 
Bn, 


The  last  step  follows  from  the  definition  of  the  derivative  which 
states  that  the  same  limit  mnst  be  obtained  for  an  approach  from  the 
right  or  the  left.  Hence  31np(nt;  a)/dnt  is  an  odd  function  as 
asserted.  Since  each  term  in  the  sum  of  (25)  has  been  shown  to  be 
zero,  the  conclusion  of  the  theorem  that  I^a(p,  a)  =  0  follows. 


Theorem  2:  Assume  that  N  samples  of  &  process  x^.  =  \is^  +  n^  are 
observed,  where  the  noise  samples  are  independent  and  each  one  is 
distributed  according  to  a  symmetric  PDF,  p(n;  a).  The  PDF  of  the 
noise  depends  upon  a  vector  a  of  unknown  parameters.  Then  the  GLRT 
(10)  for  testing  \i  =  0  vs  \i  f6  0  has  asymptotic  performance  given  by 
(15)  and  (16)  with  X  given  by  (18)  and  the  performance  is  as  good  as 
the  performance  of  the  clairvoyant  GLRT  (a  known)  as  given  by  (17). 


Proof:  follows  directly  from  previous  discussion. 
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The  next  theorem  concerns  the  implementation  of  the  GLRT  via  the 
Rao  test. 


Theorem  3:  All  statements  hold  true  in  Theorem  2  if  the  GLRT  is 
replaced  by  a  Rao  test  which  is  defined  as 


V 


N-l  . 

2  31np(  -  pst;  a) 

.  t=0  6}i 

o 

ii 

n 

(26) 


I  (0,  a) 
MM 


> 


where  a  is  the  MLE  under  Hq  and 


N-l 

I  (0,  a)  *  2  s* 


t=0 


'p'(n;  a)' 
p(n;  a)  y 


p(n;  a)dn 


(27) 


p'(  n;  a)  denotes  the  derivative  with  respect  to  n.  The  quantity  in 
brackets  is  the  Fisher  information  for  a  shift  in  mean  or  the  intr ins ic 
accuracy  of  a  PDF, ^ 


Proof:  The  asymptotic  equivalence  of  the  Rao  test  to  the  GLRT  is 
proven  in  [33 •  We  now  prove  that  the  Rao  test  for  our  problem  is  given 
by  (26)  and  (27).  Starting  with  (7)  let  §r  =  \l  and  §s  =  a  .  Also  note 
that  =  0  .  Then, 
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31np(  x;  \i,  a) 
3^i 


lr  = 


-.f 


I  (0,a) 


I  (p,  a)  is  found  as  follows. 


I  (p,  a) 


=  E 


N-l 

31n  T  p(x  -  (is  ;  a) 
t=0  *  C 

dp 


=  E 


/N-l  31np(x  -  (is  ;  a)  \2 

I  -  ' 

\  t=0 


d(i 


N-l 

=  1  E 

t=0 

N-l 

=  2  s*  E 

t=0 


31np(xt  -  pst;  a)\2' 
3(i 

f  31np(nt;  a)  \  2 
,  3n~  ) 


N-l 

=  2  S1 
t=0 

N-l 

=  I  ! 

t=0 


p ‘ (nt;  a) 

- ^  p(n  :  a)dn. 

p(nt;  a)  t  -  t 


p ' (n;  a) 
p (n;  a) 


p (n;  a)dn 


which  is  independent  of  p  and  hence  1^(0,  a)  follows, 


To  avoid  the 


integration  required  to  compute 


it 


is 
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possible  to  use  the  asymptotic  equivalence 


I  (p,  a) 

pp 


N-l 

t=0 


2 


E 


31np (n^ ; 
3n, 


2 


1 

N 


N-l 

1 

n=0 


2 

s 

n 


N-l 

2. 

t=0 


p'  (nt?  a) 
p(nt;  a) 


2 


which  for  p  =  0  (so  that  xt  =  nt)  becomes 


1(0, a) 


1 

N 


N-l 

1 

n=0 


2 

s 

u 


N-l 

I 

t=0 


P ' 


a) 

a) 


2 


(28) 


This  result  follows  from  the  law  of  large  numbers.  In  essence,  the 
integral  of  (27)  is  replaced  by  a  Monte  Carlo  evaluation.  Rao's  test 
becomes 


N-l 

I 

t=0 


91np(xt  -  ps^;  a) 

dp 


P  =  0  J 


R 


The 

21nL 


,  N-l  N-l 

j  !  i‘  i 

^  n=0  n  t=0 


p'(xt;  a) 
p(xt;  cl) 


asymptotic  statistics  of  L^'  are  identical 
q.  In  the  next  section  an  example  of  the  use 


(29) 


to  those  of  and 

of  L^'  is  given. 
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V 


Application  to  Detection  in  Non-Gauss ian  Noise 


Consider  the  problem  of  (8)  when  n^  has  a  Gaussian  mixture  PDF  or 


p  (n) 


1-8 


3 

n 


(30) 


where  <Tj  >>  ag.  The  first  term  represents  the  usual  background  noise 

while  the  second  term  is  an  interference  component.  As  an  example  a 

realization  of  the  noise  time  series  for  Og2  =  1,  Oj2  =  100  and  s  =  0.1 

is  shown  in  Figure  1.  The  noise  spikes  are  obviously  due  to  the 

contaminating  component.  The  unknown  noise  PDF  parameter  is  assumed  to 

2 

be  the  mixture  parameter  e.  It  is  assumed  that  the  noise  variances  Og 
and  Oj2  are  known.  The  known  signal  $t  is  assumed  to  be  a  DC  signal 
with  amplitude  one  so  that  the  problem  is  to  detect  a  DC  signal  \ i  of 
unknown  level  in  independent  and  identically  distributed  contaminated 
Gaussian  noise  with  an  unknown  mixture  parameter  e.  The  noise  PDF  is 
symmetric  for  all  e  in  the  range  [0,1]  so  that  the  results  of  Theorems 
1-3  apply.  The  Rao  test  of  (29)  becomes 


N-l  31np(x  -  p;  e) 

1  Jo  * 

n=o  . 

2 

(31) 

N-l 

I 

t=0 

p ' (xt;  s) 

‘  2 

p(*t;  e) 

where  e  is  the  MLE  of  e  assuming  Hq  is  true  or  e  is  the  value  of  e  over 
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the  interval  [0,1]  which  maximizes 


N-l 

J0  [  (1  -  *>  W  +  8  "l<V  ] 


(32) 


1  -  ^(n/ffg)*  -  ^(n/<r  )* 

where  0g(n)  =  - - -  e  and  0T(n)  = - -  e 


^2ji 


B 


^2n 


Some  simplifications  to  (31)  are  possible. 


dlnpti^-  ji;  e) 


3p 


H=0 


-(1  -  s) 


V  ** 


B 


(x*_  ^  ~  e  ( ~~  j?ti<V|x) 


B  t 


H  -  o 


p(*t  -  p;  e) 

[  (1  -  *)  (xt/gB}  W  +  «  Ut/qI}  W  ] 

<1  “  8)  +  8 


so  that 


4  = 


N-l  p'(x  ;  a) 

I 


.  t-0  P<V  s)  . 


N-l 

I 

t=0 


p’(xtS  e ) 
— pTx^TTT 


(33) 


where  p(xt;  a)  =  (1  -  a)  0g(xt)  +  e  0j(xt) 

p'(xt;  a)  =  -[(1  -  a)  (xt/crg2)  0g(xt)  +  a(xt/aj1)  0j(xt)] 
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The  asymptotic  performance  of  L^'  is  given  by  (15)  and  (16)  where 

\ 

X  =  H*  1^(0,  a) 

From  (27)  with  st  =  1 

f  «  f  p' (n;  e)  12 

X  =  N (i  - 7 - r  p(n;  8)dn  (34) 

J-®  L  p (n;  e)  J 

The  noncentrality  parameter  X  cannot  be  evaluated  in  closed  form  so 
that  a  numerical  evaluation  is  required. 

VI  Computer  Simulation  Results 

The  performance  of  the  detector  described  in  previous  section  is 
now  examined  via  a  computer  simulation.  To  summarize  the  assumptions 
the  detection  problem  addressed  is 

H0:  xt  =  nt 

Hl:  xt  =  v-  +  nt 

for  t  =  0,  1,  .  N-l.  N  was  chosen  to  be  500.  The  signal  level  \ i  is 
unknown  and  can  be  positive  or  negative.  The  noise  samples  nt  are 
independent  and  identically  distributed  with  PDF 
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+ 


8 


e 


•  2  2  2 
The  noise  variances  Og  and  Oj  are  assumed  to  be  known  and  =  1, 

aj  =  100  for  purposes  of  the  simulation.  The  mixture  parameter  e  is 


assumed  to  be  unknown  and  8  =  0.1.  The  SNR  is  defined  to  be 


SNR  =  10  log10  dB 

°B 


(35) 


which  represents  the  output  SNR  of  a  matched  filter.  \i  is  set  for  each 
SNR  by  using  (35).  The  Rao  test  is  given  by  (33).  A  clairvoyant 
detector  is  defined  to  be  identical  to  (33)  except  that  the  known  value 
of  s  is  used  in  place  of  its  MLE  under  Hq.  The  theoretical  asymptotic 


performance  of  the  Rao  test  and  the  clairvoyant  detector  are  given  by 


(15)  and  (16)  with  X  given  by  (34).  For  a  PpA  of  10~2  the  simulation 
results  for  the  Rao  test  and  the  clairvoyant  test  as  well  as  the 
theoretical  asymptotic  performance  are  shown  in  Figure  2.  As  expected 
the  theoretical  asymptotic  results  appear  to  accurately  predict 
performance.  The  Rao  test  performance  is  nearly  identical  to  that  of 
the  clairvoyant  test  in  accordance  with  Theorem  3. 

As  a  benchmark  to  which  performance  might  be  compared  it  is 
illustrative  to  consider  a  matched  filter.  The  matched  filter  which 
would  result  if  the  contaminating  Gaussian  PDF  were  ignored  is  to 
compare  the  sample  mean  to  a  threshold  or 


24 


The  matched  filter  requires  the  additional  knowledge  of  the  sign  of  (i. 
For  the  simulation  p  was  chosen  to  be  positive.  The  matched  filter 
might  be  used  if  the  mixture  parameter  were  unknown.  This-  would 
correspond  to  assuming  that  e  =  0  when  in  fact  s  >  0.  A  limiter  of 
course  would  improve  the  performance.  The  results  are  shown  in  Figure 
3.  The  Rao  test  allows  one  to  detect  signals  30  dB  lower  than  a 
matched  filter,  clearly  illustrating  the  importance  of  using  the 
mixture  parameter  information  in  any  detector. 

VII  Conclusions 


We  have  shown  that  for  detection  of  known  signals  of  unknown 
amplitude  in  non-Gaussian  noise  whose  PDF  is  not  completely  specified 
that  near  optimal  performance  may  be  obtained.  The  optimality  is  an 
asymptotic  optimality  (large  data  records  or  small  signals)  but 

generally  holds  for  problems  of  practical  interest.  The  optimal 

t 

detector  is  obtained  by  using  a  GLRT  or  a  Rao  test,  the  latter  being 
easier  to  implement,  and  is  based  on  the  assumption  that  the  PDF  of  the 
noise  is  symmetric.  This  condition  is  not  overly  restrictive  in  that 
nearly  all  real  world  noise  has  this  property.  Another  benefit  of  the 
theory  is  that  the  detection  performance  is  easily  computed 
analytically.  Future  work  will  extend  the  results  to  the  case  of 
non-independent  noise. 
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Figure  2  -  Comparison  of  Detection  Performance  - 
Rao  Test'  and  Clairvoyant  Rao  Test 
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Figure  3  —  Comparison  of  Detection  Performance  — 
Rao  Test  and  Matched  Filter 
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Appendix  A 


Example  of  GLRT  Clairvoyant  Performance 

In  this  appendix  the  GLRT  for  two  detection  problems  are  derived. 
Consider  the  generic  problem  of 

H0  :  xt  =  nt 
Hi  :  xt  =  (i  +  nt 
for  t  =  0,  1,  . .  . ,  N-l. 


is  an  unknown  amplitude  either  positive  or  negative.  nt  is 
independent  and  identically  distributed  noise  samples  with  each  sample 
having  the  PDF  ~  N(0,  a2).  In  problem  1  the  noise  variance  a2  is 
assumed  unknown  while  in  problem  2,  cr2  is  known  a  priori.  The  GLRTs 
for  problems  1  and  2  respectively  are 
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N-l 
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t=0  ^77 


It  is  well  known  that  for  problem  1 
N-l 
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t=0 

N-l 
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t=0 
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A  1 
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a3  1 

°  =  N 


X^  =  X 


I  (xt  -  x)‘ 
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AS  1  1 

«  =  «  I 
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t=0 


and  that  for  problem  2 


(A-2) 


Substitution  of  the  MLEs  into  (A-l)  and  (A-2)  yields 


l„  -  (j'/;')N/2 


lgc 


2a 


N-l  N-l 

I  (x  -  x)*  -  I  x* 
t=0  t=0 


0  .  a2  *2  _ 2 
Since  a  =  a  -  x  , 
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cr 


so  that  finally 
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Also,  after  some  simplification 
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which  yields 


2  ln  L, 


OFFICE  OF  NAVAL  RESEARCH 
STATISTICS  AND  PROBABILITY  PROGRAM 

BASIC  DISTRIBUTION  LIST 
FOR 

UNCLASSIFIED  TECHNICAL  REPORTS 


FEBRUARY  1982 


Copies 


Statistics  and  Probability 
Program  (Code  411 (SP)) 

Office  of  Naval  Research 
Arlington,  VA  22217  3 

Defense  Technical  Information 
Center 

Cameron  Station 

Alexandria,  VA  22314  12 

Commanding  Officer 
Office  of  Naval  Research 
Eastern/Central  Regional  Office 
Attn:  Director  for  Science 
Barnes  Building 
495  Summer  Street 
Boston,  MA  02210  1 

Commanding  Officer 
Office  of  Naval  Research 
Western  Regional  Office 
Attn:  Dr.  Richard  Lau 
1030  Eas-t  Green  Street 
Pasadena,  CA  91101  1 

U.  S.  ONR  Liaison  Office  -  Far  East 
Attn:  Scientific  Director 
APO  San  Francisco  96503  1 

Applied  Mathematics  Laboratory 
David  Taylor  Naval  Ship  Research 
and  Development  Center 
Attn:  Mr.  G.  H.  Gleissner 
Bethesda,  Maryland  20084  1 

Commandant  of  the  Marine  Coprs 
(Code  AX) 

Attn:  Dr.  A.  L.  Slafkosky 
Scienti fic  Advisor 
Washington,  DC  20380  1 


Copies 


Navy  Library 

National  Space  Technology  Laboratory 

Attn:  Navy  Librarian 

Bay  St.  Louis,  MS  39522  1 

U.  S.  Army  Research  Office 
P.0.  Box  12211 
Attn:  Dr.  J.  Chandra 
Research  Triangle  Park,  NC 

27706  1 

Di rector 

National  Security  Agency 

Attn:  R51 ,  Dr.  Maar 

Fort  Meade,  MD  20755  1 

ATAA-SL,  Library 
U.S.  Army  TRADOC  Systems 
Analysis  Activity 
Department  of  the  Army 
White  Sands  Missile  Range,  NM 
88002  1 

ARI  Field  Unit-USAREUR 

Attn:  Library 

c/o  ODCSPER 

HQ  USAEREUR  &  7th  Army 

APO  New  York  09403  1 

Library,  Code  1424 
Naval  Postgraduate  School 
Monterey,  CA  93940  •  1 

Technical  Information  Division 
Naval  Research  Laboratory 
Washington,  DC  20375  1 

OASD  { I&L) ,  Pentagon 
Attn:  Mr.  Charles  S.  Smith 
Washington,  DC  20301 


1 


Copies  Copi 


Di rector 
AMSAA 

Attn:  DRXSY-MP ,  H.  Cohen 
Aberdeen  Proving  Ground,  MD  1 
21005 

Dr.  Gerhard  Heiche 
Naval  Air  Systems  Command 
(NAIR  03) 

Jefferson  Plaza  No.  1 
Arlington,  VA  20360  1 

Dr.  Barbara  Bailar 
Associate  Director,  Statistical 
Standards 
Bureau  of  Census 
Washington,  DC  20233  1 

Leon  Slavin 

Naval  Sea  Systems  Command 
(NSEA  05H) 

Crystal  Mall  M,  Rm.  129 
Washington,  DC  20036  1 

B.  E.  Clark 

RR  § 2,  Box  647-8 

Graham,  NC  27253  1 

Naval  Underwater  Systems  Center 
Attn:  Dr.  Derrill  J.  Bordelon 
Code  601 

Newport,  Rhode  Island  02840  1 

Naval  Coastal  Systems  Center 
Code  741  . 

Attn:  Mr.  C.  M.  Bennett 
Panama  City,  FL  32401  1 

Naval  Electronic  Systems  Command 
(NELEX  612) 

Attn:  John  Schuster 
National  Center  No.  1 
Arlington,  VA  20360  1 

Defense  Logistics  Studies 
Information  Exchange 
Army  Logistics  Management  Center 
Attn:  Mr.  J.  Dowling 
Fort  Lee,  VA  23801  1 


Reliability  Analysis  Center  (RAC) 

RADC/RBRAC 

Attn:  I.  L.  Krulac 

Data  Coordinator/ 

Government  Programs 
Griffiss  AFB,  New  York  13441 

Technical  Library 
Naval  Ordnance  Station 
Indian  Head,  MD  20640 

Library 

Naval  Ocean  Systems  Center 
San  Diego,  CA  92152 

Technical  Library 
Bureau  of  Naval  Personnel 
Department  of  the  Navy 
Washington,  DC  20370 

Mr.  Dan  Leonard 
Code  8105 

Naval  Ocean  Systems  Center 
San  Diego,  CA  92152 

Dr.  Alan  F.  Petty 
Code  7930 

Naval  Research  Laboratory 
Washington,  DC  20375 

I 

Dr.  M.  J.  Fischer 
Defense  Communications  Agency 
Defense  Communications  Engineering 
Center 

1360  Wiehle  Avenue 
Reston,  VA  22090 

Mr.  Jim  Gates 
Code  9211 

Fleet  Material  Support  Office 
U.  S.  Navy  Supply  Center 
Mechanicsburg,  PA  17055 

Mr.  Ted  Tupper 
Code  M- 31 1 C 

Military  Sealift  Command 
Department  of  the  Navy 
Washington,  DC  20390 


Copies 


Copies 


Mr.  F.  R.  Del  Priori 
Code  224 

Operational  Test  and  Evaluation 
Force  (OPTEVFOR) 

Norfolk,  VA  23511  1 


